AD-A179  444 


REPORT  DOCUMENTATION  PAGE 


Unc lass  i  f ied 


W  r  .  V.M  t  HJIN 


12a.  SECLR'Tv  CLASSIFICATION  A  JMi 


2b.  0  £  C  L  ASS  I  F  i  C  A  T  i  ON  OOWNGRA1 


4  PERFORMING  ORGANIZATION 


6a.  NAME  OF  PERFORMING  ORGANIZATION  |6o.  OFFICE  SYMBC 

The  University  of  Maryland  'ir applicable) 


lb.  RESTRICTIVE  MARKINGS 


|  3.  DiSTRlBUTION/AVAl  LABILITY  OF  REPORT 

Approved!  for  public  ?»Xm997 
distribution  unlimited 


5.  MONITORING  ORGANIZATION  REPORT  NUMBERlSI 

AFOSR  TW-  87-  0  477 


7a.  NAME  OF  MONITORING  ORGANIZATION 

AFOSR 


6c.  address /City,  state  and  /ip  co<u>  Oniv,  of  Maryland 
College  Park  Campus 

Laboratory  for  Plasma  and  Fusion  Energy  Stuc 
Energy  Research  31dg,  College  Park  Md  20742 


•«.  NAME  OF  FUNOING/SPONSORING 
ORGANIZATION 

AFOSR 


8c  AOORESS  City.  State  and  ZIP  Code) 


same  as  7b . 


11  tiTlE  inciude  Secunty  Classification/ 


8b.  OFFICE  SYMBOL 
•If  applicable  > 


7b.  AOORESS  (City,  State  and  ZIP  Code) 

ies/  Bld§  4*0 

Bolling  AFB ,  D.C.  20332-6448 


9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 

AFOSR-86-0055 


10  SOURCE  OF  FUNDING  NOS. 


WORK  UNIT 
NO 


"FORMATION  AND  STABILITY  OF  PARTIALLY-NEUTRALIZED  PLASMA  CLOMPS" 


12  PERSONAL  AUTMORlSl 

Dr  .  John  Gui 1 lorv 


13a  type  OF  REPORT 

Annua  1 


130.  time  covereo 

gHO^6/01/01  r086/ 12/01 


14  OATE  OF  REPORT  lYr  .  .Mo..  Day) 

87  January  1 


15  page  count 

72 


COSATI  COOES 


18  SUBJECT  TERMS  'Continue  on  r^ier te  if  necessary  and  identify  by  block  number t 

collective  Accleration,  Plasmoids'  Electron  Beams 


19  ABSTRACT  Continue  on  reverse  if  necessary  and  identify  by  block  number* 

Experimentally  a  University  of  Maryland  group  has  investigated  neutralized  electron 
beam  propagatio  in  nonconducting  vacuum  drift  tubes  (to  eliminate  image  currents  and  a 
clear  return  path  for  the  electrons)  and  investigate^!  a  significant  new  way  of  in¬ 
creasing  the  ion  velocity  (and  therefore  the  propagation  distance)  in  such  systems. 

In  thisprogram  the  investigators  have  constructed  a  microscopic  equilibrium  for  such  a 
propagating  electron/ion  beam  in  which  both  species  are  in  radial  force  balance,  and 
they  have  begun  numerical  simulation  studies  of  bean  and  ion  propagation  in  such  a  system. 


20  OiSTRi  BUTlON.  A  V  Al  LABI  LI  TY  OF  ABSTRACT 


UN  CL  ASSI  F  >  E  O  UNLIMITEO  _  SAME  AS  RPT  CC  OTiC  USERS  _ 


22<  NAME  OF  RESPONSIBLE  iNDiviOual 


BRUCE  L.  SMITH 


21  ABSTRACT  SECURITY  CLASSIFICATION 

UNCLASSIFIED 


22b  TELEPHONE  NUMBER 
Include  \reaCodei 

(202)767-^907 


22c  OFFICE  SYMBOL 


m 


afosr 


T»-  8  7-0  477 


ANNUAL  TECHNICAL  REPORT 

FORMATION  AND  STABILITY  OF  PARTIALLY-NEUTRALIZED 
PLASMA  CLUMPS 

Grant  No.  AFOSR-86-0055 

For  the  period  -‘Throve  ’  ^  •  - 

January  1  ,  1986  to  December  15,  198flistri  >  • 

j.o:j 

Submitted  to 

Air  Force  Office  of  Scientific  Research 
January  1987 


Prepared  by 

Laboratory  for  Plasma  and  Fusion  Energy  Studies 
University  of  Maryland 
College  Park,  MD  20742 


/».■;  'rOPr.c 

I 

*  A  i  t  ji. 

-,V;  r«l  - . '  T  ,• 

M.V  Tt.-fV.. 


V:  •  *•  ’  '"[T  •v-JCAflCH  (ATSC) 
1  ’  '■  ■  ;v. 'awed  and  is 

"i;cn 


•*  i.  v' 


RESEARCH  PROGRESS  AND  FORECAST  REPORT 


Submitted  to: 
Submitted  by: 


Grant  Number: 


Air  Force  Office  of  Scientific  Research 

Laboratory  for  Plasma  and  Fusion  Energy  Studies 
University  of  Maryland 
College  Park,  MD  20742 

AFOS  R-86-0055 


Period : 


January  1,  1986  to  December  15,  1986 


Principal  Investigators:  John  Guillory,  Senior  Research  Associate 

Martin  P.  Reiser,  Professor 
Charles  D.  Striffler,  Professor 

Title  of  Research  Project:  Formation  and  Stability  of  Partially- 

Neutralized  Plasma  Clumps 


Acce-jion  lor 

— r-> 

NTIS  CR4&I 
D1IC  TAB 

Unannounced 

Justification 

U 

u 

by 

U.it  ibulic/ 

~~~  ~  '  -  ~  - - -i 

. .  . i 

Avoilabilit 

y  Ce  lts 

Dist 

Avji1  _  •'! :  or  ; 

SpvCiJ! 

O'! 

.  —  ■  - - 

jf/L  *-<.  -S  *  a  -Nil  ».  k  /  j, , 


PROJECT  SUMMARY 


PROJECT  TITLE:  Formation  and  Stability  of  Partially-Neutralized  Plasma 
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Theoretical  and  numerical  studies  of  plasmoid  shedding  from  gas  puffs  and 
virtual  cathodes;  and  stability  of  plasma  clumps  traversing  low-density 
plasma  background. 

APPROACH: 

Analytic  and  numerical  analysis  of  time-dependent  electron  and  ion 
emission  from  downstream  side  of  a  plasma  puff  subjected  to  an  electron 
beam  pulse;  analytic  and  numerical  analysis  of  equilibrium  and  stability 
of  a  bounded,  relatively  dense  plasma  in  a  tenuous  background  by  treating 
the  background  as  "beam"  in  the  plasma  frame. 

This  report  describes  the  first  year  of  the  suggested  two-year  effort. 
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1.  INTRODUCTION 


Studies  of  the  propagation  of  Intense  charged  particle  beams  into 
vacuum  have  been  pursued  at  the  University  of  Maryland  under  AFOSR 
sponsorship  since  1984.  That  program  has  been  investigating  under  what 
conditions  intense  charged  particle  beam  energy  may  be  propagated  into 
free  space  without  the  benefit  of  confining  magnetic  fields  or  conducting 
boundaries.  Experimental  and  theoretical  work  has  centered  around  the 
propagation  of  intense  electron  beams  into  vacuum  after  passage  through 
an  ion  source  localized  near  the  injection  point.  In  such  systems 
propagation  can  occur  when  ions  are  accelerated  from  the  source  region 
into  vacuum  by  the  space  charge  fields  of  the  electron  beam.  The 
resultant  neutralization  of  the  electron  beam  by  the  accelerated  ions 
allows  effective  propagation  in  the  vacuum  region,  to  the  extent  studied 
in  the  experiments.  The  goal  of  that  related  work  is  the  generation  of 
nearly  charge-neutral,  current-neutral  charged  particle  beams  which  are 
capable  of  propagating  beam  energy  in  vacuum  at  a  significant  fraction  of 
the  velocity  of  light. 

From  the  related  work  in  our  group,  we  have  arrived  at  the  following 
understanding  of  the  dynamics  when  an  intense  electron  beam  is  injected 
into  vacuum  through  a  localized  plasma  source: 

(a)  At  the  beginning  of  electron  beam  injection,  a  virtual  cathode 
is  established  in  the  region  of  the  localized  plasma  source.  The 
electric  field  associated  with  the  virtual  cathode  reflects  the  electrons 
and  accelerates  the  ions  to  a  mean  kinetic  energy  equal  to  the  electron 
beam  energy,  cVq. 


V. 


(b)  After  acceleration,  the  bulk,  of  the  ions  drift  with  a  mean 
velocity  v^  =  (2eVg/mi)^“. 

(c)  As  a  quasi-steady  ion  flow  is  established,  the  virtual  cathode 
moves  a  short  distance  downstream  and  then  remains  at  this  new  position. 

(d)  During  the  motion  of  the  virtual  cathode  to  its  new  position,  a 
group  of  ions  in  phase  with  the  moving  potential  well  gains  additional 
kinetic  energy  qEAz  where  E  is  the  average  electric  field  associated  with 
the  virtual  cathode  and  Az  is  the  displacement  distance  of  the  virtual 
cathode.  This  phenomenon  represents  an  explanation  for  the  so-called 
"high-energy  tail"  in  the  ion  distribution  observed  in  the  experiment. 

(e)  A  large  fraction  of  the  injected  electrons  are  reflected  at  the 
displaced  virtual  cathode.  The  rest  propagate  together  with  the  ions 
that  form  a  charge-neutralizing  channel.  The  equilibrium  state  of  the 
self-focused  electron-ion  beam  propagating  downstream  has  been  studied. 
The  initial  results  of  that  study  are  presented  in  a  paper  to  be 
published  in  the  proceedings  of  the  Beams  86  Conference  (see  Appendix 
A).  Improved  results  were  presented  at  the  1986  Plasma  Physics  Meeting 
(see  Appendix  A). 

(f)  It  appears  that  a  second  virtual  cathode  forms  at  the  beam  front 
from  which  most  of  the  propagating  electrons  are  reflected. 

The  purpose  of  this  related  study,  begun  in  January  1986,  is  to 
model  and  analyze  the  formation  of  partially  neutralized  plasma  clumps 
emitted  from  such  an  ion-source  plasma  traversed  by  an  electron  beam,  and 


3 


2.  FORMATION  -  ID 


2.1  Ion  Motion 

Previous  1-D  particle  codes  [Dantsker  et  al..  Bull.  APS  29_,  1353 
(1984)]  have  indicated  that  most  of  the  ions  (and  electrons)  follow  each 
other  along  a  single  phase-space  trajectory.  This  implies  that  a  fluid 
model  may  be  adequate  for  a  1-D  description  of  the  accelerated  motion  of 
the  plasma  ions  as  they  leave  the  ionizing  source  plasma.  To  this  end, 
we  then  begin  with  the  following  simple  1-D  fluid  model  for  the  ions: 


3  n  +3  (n  v  )  =  S  (z )  , 
t  1  z  1  i  i 


3  zE  =  4Tre(n,  -  n^)  , 

-z/i 

with  n^(t  =  0)  =  n^e  >  £  “  gas  puff  scale-length.  Here  S^(z)  is  the 

ionization  source  rate,  and  n^  must  be  coupled  to  E. 

To  simplify  the  expansion  of  this  ion  plasma  even  further,  let  us 
assume  £  is  small,  neglect  the  production  of  new  ions  (i.e.,  take  = 
0)  ,  and  look  for  a  self-similar  expansion.  Letting  x  =  z/vgt  be  the 
similarity  variable,  where  vq  is  the  ion  velocity  when  given  the  input 
beam-electron  energy  E^.  Rewriting  the  above  equations  in  the 


n  =  n^/  (n^  at  4)  =  0)  , 


and  using  the  similarity  condition,  one  arrives  at 


.  ,  dn  du 

(u  -  x)  -7—  +  n  -j— 
dx  dx 


=  0 


d'l> 

dx 


+ 


(u  -  x)  ~  =  0 
dx 


and  we  require  some  form  n(4<)  for  the  ions  and  beam  (the  resulting  plasma 
expansion  is  quasineutral).  For  simplicity  of  results  we  assume 


n(40  =  (1  +  40P 


with  p  an  adjustable  power  (typically,  experiments  indicate  p  ~  3).  The 
solution  of  this  simplified  set  of  ODE's  is 


u  =  Uq  +  ax  , 
n  =  pp[uQ  -  (1  -  a)x]2p 
'P  =  p[uQ  -  ( 1  -  a  )x] 2  -  1  , 

with  a  =  1  -  (2p  +  l)-1/2  and  u0  =  21,/2(2p  +  l)-1/2,  from  ion  energy 
conservation,  u2/2  +  4*  *■  const.  =  0  at  x  =  0  out  to  x  =  u^/(l  -  a)  =  *  7, 


where  n  goes  to  zero. 


The  peak  Ion  kinetic  energy,  1/2  u^  at  x  =  Uq/(1  -  a)  ,  Is  unity, 
i.e.,  ions  at  the  head  of  the  expansion  have  the  full  electron  beam-input 
energy.  When  electron  energy  loss  to  the  potential  is  included,  the 
value  of  a  increases  and  the  peak  ion  kinetic  energy  can  be  somewhat 
larger. 

The  transient  introduction  of  the  electron  beam  leads  to  an 
oscillatory  transient,  damping  to  such  a  self-similar  expansion.  Such 
transients,  the  subject  of  part  of  our  second-year  effort,  are  intimately 
related  to  the  electron  dynamics,  and  are  Important  if  large,  because 
they  may  determine  the  size  of  emitted  plasmoids. 


2.2  Electron  Motion 


To  model  the  beam  and  sheath  electron  disturbances  at  the  exit- 
surface  of  a  relatively  cold  dense  plasma  traversed  by  the  electron  beam, 
we  begin  with  the  1-D  model  system 


9 2 (ed> )  =  ATTe2(n  +  n  -  n  )  , 
z  e  b  1 


9  n  +  9  (n  v  )  =  S  , 
t  e  zee  e 


,  ej> , 


9v  +  v9v  =  3  ( — ) 
t  e  e  z  e  z  m 


Here  we  let  T  =  oj^OH,  x  =  z/hn  where  hn  is  the  "initial"  or  typical 


density  gradient  scale  height,  and  u  =  v>^1!Jpy5Q•  berms  of  these 


dimensionless  variables  and 
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the  above  equations  become 


V  +  I(*)  , 


3tu 


3  >  , 
x 


3  N'  +  u3  N  +  N’  3  u  =  S  , 
i  e  x  e  ex 


where  we  have  neglected  v^  •  Vve  compared  with  3tve. 

This  set  of  equations  purports  to  describe  space  charge  oscillations 
in  the  non-quasineutral  sheath  of  the  quasi-stat lonary-ion  plasma,  i.e., 
it  is  intended  as  a  fast-timescale  fluid  description  of  the  electron 
sloshing.  The  function  I(>>)  simply  describes  the  profile  in  x  of  the 
beam  electron  density,  and  is  normalized  to  unity  at  x  =  0.  For  a 
uniform  beam  (reflexing  or  not),  1=1.  ^  is  specified,  in  a  way 
consistent  with  the  ion  slow-timescale  equations,  e.g.,  the  self-similar 
model  treated  above,  or  some  more  general  one. 

Analysis  (and/or  generalization)  of  this  1-D  electron-motion/space- 
charge  oscillation  model  is  part  of  our  projected  second-year  effort,  as 
is  computer  simulation.  When  the  electron  beam  is  modulated  in  time 
(even  if  only  by  its  rise  and  decay)  we  expect  that  the  large-amplitude 
space-charge  oscillations  near  the  forward  edge  of  the  quasineutral 


expanding  plasma  can  lead  to  modulation  of  the  ion  flow.  We  cannot  yet 


predict  the  relative  phases  of  the  ion  and  electron  modulation 
(jetting?),  although  energetics  would  tend  to  favor  in-phase 
(quasineutral)  clump  formation. 

In  Fig.  1  a  schematic  of  the  model  described  in  Secs.  2.1  and  2.2  is 
shown.  The  electron  beam  is  shown  incident  onto  the  plasma  source  region 
and  the  resulting  virtual  cathode  forming  just  beyond  the  downstream 
plasma  surface.  An  early  time  profile  of  the  potential  is  displayed  and 
the  indicated  movement  of  the  potential  profile  to  the  downstream  plasma 
surface.  The  sheath  region,  from  which  ions  are  emitted,  is  also 
displayed  on  the  downstream  plasma  surface.  Ions  are  drawn  out  of  this 
region  and  form  schematically  the  indicated  spatial  profile. 


2.3  1-D  Simulations 


The  numerical  model  of  our  localized  plasma  source-beam  propagation 

system  is  shown  in  Fig.  2a.  A  solid  electron  beam  of  radius  Rb,  current 
2  2 

Iq  =  and  voltage  eVp  =  me  (Yq  -  1)  is  injected  through  one  end  of 

an  evacuated  drift  tube  of  radius  Rw  and  length  L.  The  entire  system  is 
immersed  in  a  magnetic  field  that  is  sufficiently  large  to  impede  any 
radial  expansion.  Without  the  presence  of  a  source  of  ions,  a  deep 
potential  well  of  depth  about  equal  to  the  voltage  Vq  forms  in  a  region 
about  1-2  mm  from  the  entrance  end  when  the  injected  current  is  above  the 
space-charge  limiting  current 

2/3  3/2 

17(Yq  -  ir' 

XL  =  1  +  2in(R  ART)  "  kA  ' 

w  b 


In  this  case  the  majority  of  the  incoming  beam  particles  are  reflected 


Sheath  region: 
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back  to  the  anode  (z=0)  by  this  deep  potential  well  and  a  fraction  are 
j  transmitted  past  it. 

In  Fig.  2b,  we  have  pictorially  shown  where  plasma  is  formed  due  to 
the  presence  of  a  cloud  of  gas  localized  to  the  region  near  the  entrance 
end.  This  cloud  of  gas  is  ionized  by  the  electron  beam  and  subsequently 
by  the  plasma  particles  (ions  and  electrons)  and  serves  as  the  source  of 
ions  that  enables  the  beam  to  propagate  further  downstream.  The 
|  development  of  the  code  (J.  Grossmann's  Ph.D.  thesis,  1982)  and  a 

prescribed  rate  of  plasma  formation  had  been  the  model  studied  during  the 
initial  phases  of  research.  In  1984,  we  quantified  the  ionization 

process.  This  model  is  described  below  and  results  shown  of  the  initial 
runs  with  our  electrostatic  code  [(r,z)  geometry;  pz  momentum]. 

A  neutral  gas  is  assumed  to  occupy  a  region  close  to  the  anode. 
Specifically 

I 


;  P(z)  = 

! 


0  <  z  <  Zq 


z0  <  2  <  L 


> 


where  p(z)  is  the  pressure  of  the  gas  in  torr.  The  formation  of  plasma 
particles  via  ionization  is  given  in  terms  of  a  cross-section  of  the 
interaction.  Mathematically  the  formation  of  plasma  particles  (ions  and 
electrons)  in  a  time  At  is  given  by 


an  =  an  =  l 
pl  Pe  a 


Sa ( KEa )  X  va  X 


At  x 


Na  particles 


where  Sa  is  the  number  of  ionization  events  per  cm  per  torr,  v^  and  KEa 
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are  velocity  and  kinetic  energy  of  the  ionizing  particle  resulting  in  the 
ionization,  and  is  the  number  of  particles  of  type  a.  The  sum  over 
species  a  represents  the  contribution  to  ionization  by  beam  electrons, 
accelerated  plasma  ions,  and  plasma  electrons.  The  time  At  is  chosen  as 
the  time  step  of  the  simulation  and  AN  is  calculated  for  each  mesh 
region  Az  separately.  When  a  sufficient  number  of  plasma  particles  have 
been  formed  in  a  mesh  region  a  macroparticle  of  that  type  is  released 
into  the  system  and  is  subsequently  accelerated  by  the  electric  field. 
This  particle  is  then  also  allowed  to  ionize  the  gas.  The  specific 
expressions  for  incorporated  into  the  numerical  model  are: 


KE  <  20 

e  _ _  o 

20  <  KE  <  /10  x  10 
e  - 

/TO  x  102  <  KE  <  105 

5  e 

10  <  KE 

e 


and 


KEi  <  10 

10  <  KE  0.3*  104 
A. 3  x  104  <  KE 


e 


where  KE^  Is  expressed  in  terras  of  eV  and  these  expressions  model  the 
curves  given  in  the  reference  C.  L.  Olson,  V.  Schumacher,  Collective  Ion 
Acceleration,  Springer  Tracts  #84,  1979,  p.  50. 

The  results  of  our  numerical  simulation  incorporat ing  the  ionization 
process  described  above  are  shown  in  the  set  of  plots  in  Figs.  3(a)- 
3(c).  The  system  and  beam  parameters  are  listed  in  the  figure  caption 
for  Fig.  3(a).  Note  the  following  points: 

-  The  t  =  0.3  ns  plots  show  the  potential  depression  region  due  to 

the  high  current  electron  beam.  That  is,  near  z  ~  1-2  mm  the  beam 

electrons  are  reflected  back  to  the  anode  and  only  a  few  are  transmitted 

downstream.  This  latter  number  represents  the  limiting  current 

I  ~  2  kA.  These  plots  also  show  how  and  where  plasma  particles  are 
L 

formed.  That  is,  they  are  forming  near  the  potential  minimum  since  the 
rate  of  ionization  is  largest  when  the  beam  electrons  have  a  relatively 
low  velocity. 

-  The  t  >  3  ns  plots  show  the  results  for  the  entire  tube  length. 
The  potential  depression  $  ~  -  is  restricted  to  the  anode  region  early 
in  time  but  broadens  as  time  increases.  At  the  same  time  ions  move 
downstream,  the  region  of  broadening  being  about  equal  to  the  extent  of 
the  region  occupied  by  the  ions. 

-  The  ions  at  the  front  of  the  broadened  region  have  an  oner g.'  it 
about  2  MeV,  i.e.,  four  times  the  beam  energy  in  the  diode.  The  bulk  of 
the  ions  have  an  energy  between  0.5  and  1.0  MeV.  This  increase  above  o.S 
MeV  appears  to  be  due  to  the  positive  potential  region  restricted  to  the 
gas  cloud  region.  This  positive  potential  oscillates  about  0.7  MV. 


12 


V*  * 

i 

kv. 


-  For  t  ^  4.5  ns  Che  steep  gradient  in  potential  remains  near  the 
end  of  the  gas  cloud  region  for  this  gas  pressure.  However,  though  there 
is  no  further  increase  in  the  peak  ion  energies,  the  ions  are  still 
supplied  to  the  beam  front  for  beam  propagation  purposes. 

-  By  keeping  count  of  the  number  of  plasma  particles  formed,  we  find 
that  no  movement  of  the  deep  potential  region  occurs  until  ionization  by 
the  plasma  ions  becomes  important,  i.e.,  until  avalanching  occurs. 

The  results  shown  here  represent  a  great  improvement  on  our  previous 
numerical  studies,  mainly  because  of  the  quant i f icat ion  of  the  gas  cloud 
region.  Further  studies  will  be  performed  during  the  second  year  of  this 
research. 
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3.  FORMATION  -  2D  EFFECTS 


3.1  Analytic  Theory 


A  current-carrying,  ablating  planar  plasma  is  subject  to  a  thermal- 
resistive  instability  when  inelastic-collision  smoothing  of  temperature 
gradients  is  not  dominant  (Tripathi,  Ottinger,  and  Guillory,  1983). 
Hotter  regions  carry  more  current  and  heat  further,  under  a  range  of 
circumstances.  This  leads  to  radial  (or  otherwise  transverse) 
inhoraogeneities  of  the  exit  surface  of  the  source  plasma,  on  a  timescale 
that  can  be  approximately  the  same  as  the  beam  risetime  (~  10  ns  in  the 
Maryland  experiments).  In  the  first  year  of  this  grant  research,  we  have 
coupled  the  first-order  (linearized)  perturbations  in  density  and 
temperature  to  the  axial  ion  motion,  i.e.  ,  we  have  analyzed  the  ion 
motion  response  as  driven  by  the  unstable  electron  currents.  And  we  have 
generalized  the  perturbation  equations  of  the  instability  analysis  to 
allow  for  axial  gradients  in  the  equilibrium  density,  temperature,  etc. 

In  order  to  explore  the  f ilamentation  and  jetting  of  a  plasma  "slab" 
through  which  an  electron  beam  is  passing,  we  have  derived  linear 
perturbation  equations  for  the  low-frequency  oscillations  in  a 
quasineutral  plasma  with  gradients  in  the  axial  direction.  The  rationale 
for  this  is  that  transverse  f ilamentation  and  axial  jet  structures  may 
determine  the  size  of  quasineutral  plasma  clumps  shed  from  the  downstream 
edge  of  the  plasma.  There  is  both  theoretical  and  experimental  evidence 
that  such  transverse  disturbances  may  occur  in  the  cathode  blowoff  plasma 
of  energetic  electron  beam  diodes  (Tripathi,  Ottinger,  and  Guillory, 
1983;  Stephanakis,  Goldstein,  and  Mosher,  1982). 


The  perturbation  equations  comprise  a  set  of  coupled  ordinary 
differential  equations  in  z,  with  complex  eigenf requency  w,  for  the 
perturbation  eigenfunctions  of  density,  current,  temperature,  and 
fields.  As  a  first  step  we  have  analyzed  the  ion  response  as  driven  by 
the  unstable  electron  current  perturbations. 

The  model  is  as  follows:  a  quasineutral  plasma  initially  uniform  in 

x  and  y  has  a  density  falloff  in  z  with  scale  length  k  The  z-directed 

fluid  velocity  v  in  the  plasma  has  a  scale  length  k  ^  in  z.  The  z- 

directed  ion  fluid  velocity  v^  has  the  same  scale  length,  .  There  may 

also  be  electron  temperature  gradients  in  z,  with  scale  length  k,^.  The 

ions  are  cold.  Through  this  plasma  streams  an  energetic  electron  beam, 

which  travels  with  only  (  +  and  -)  z-dlrected  motion  and  extends  uniformly 

in  z  beyond  the  plasma.  There  is  a  magnetic  field  in  the  y-direction  due 

to  the  electron  current.  We  perturb  the  equations  of  this  plasma  with  a 

perturbation  varying  as  exp(ik  x  +  ik  y  -  iwt),  but  with  phase  velocity 

x  y 

<<  c.  We  assume  the  perturbed  ion  and  electron  densities  are  equal 
(quasineutral  perturbations).  The  beam  electrons  are  assumed  unperturbed 
in  their  brief  transit  through  the  perturbations  in  the  plasma. 
Collisions,  including  inelastic  collisions  are  included  in  the  finite- 
resistivity  plasma  model,  although  they  become  negligible  where  the 
density  becomes  small.  This  is  because  we  know  from  our  previous  work 

with  a  uniform  plasma  slab  model  that  inelastic  electron  collisions 

affect  the  stability  of  the  system  (Tripathi,  Ottinger,  and  Guillory, 
1983).  The  plasma  is  assumed  to  have  a  small  charge  imbalance 
compensating  the  electron  beam  space  charge,  and  this  charge  imbalance  is 
not  negligible  where  the  density  is  low.  It  couples  with  the  electric 
field  g  to  give  a  self-consistent  force  in  addition  to  the  J  *  B  and  vp 
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forces. 

The  perturbed  quantities  are  thus  n,  p  =  -  n  ,  T(=  T0),  v, 

J(=  Jy),  E  and  B.  But  since  perturbations  of  the  form  and  give 

rise  to  high  frequency  electron  plasma  waves  to  which  the  ions  cannot 
respond  except  as  a  nonlinear  ponderomot ive  matter,  we  ignore  them  in 
this  particular  analysis.  (We  consider  ponderomot ive  ion  acceleration 
elsewhere. ) 

In  its  simplest  form,  the  thermal-resistive  instability  displays  the 
growth  rate 
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with  the  definitions: 


=  electron  collision  rate 
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inelastic-collision  energy  loss  rate  of  electrons 
(proportional  to  radiation  per  electron) 


Coupling  the  perturbations  of  F.  field  and  electron  temperature  to 
perturbed  ion  motion,  we  get,  e.g.,  for  the  perturbation  amplitude  of 
axial  ion  velocities 
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with  the  detinitions 


ionization 

and  with  E^  and  the  electric  field  and  electron  temperature 

perturbat ions . 
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3.2  2-D  Code 


We  have  loaded  onto  the  local  computer  the  2-D,  fully- 
electromagnetic  particle  code  WAVE,  originated  at  Los  Alamos  for  the 
study  of  hot  electron  (etc.)  dynamics  near  the  absorption  of  laser  light 
in  an  ablating  plasma. 

The  WAVE  code  solves  Maxwell's  equations  and  particle  equations  of 
motion  on  a  Cartesian  mesh  with  a  variety  of  field  and  particle  boundary 
conditions.  Although  it  has  been  used  primarily  to  study  the  interaction 
of  intense  laser  light  with  plasmas,  it  has  the  capability  of  studying 
many  other  phenomena  including  shock  waves,  beam  plasma  instabilities, 
etc.  It  was  first  written  by  R.  L.  Morse  and  C.  W.  Nielson  in  1971, 
improved  by  C.  W.  Nielson  and  E.  L.  Lindman,  maintained  by  D.  W.  Forslund 
(with  the  addition  of  highly  optimized  disk  I/O  routines),  and  recently 
modified  by  C.  H.  Aldrich  to  include  particle  beam  injection  and  by  D.  W, 
Forslund  to  include  aperiodic  boundary  conditions  in  the  y-direction,  to 
conform  closely  to  the  ANSI  FORTRAN  77  standard,  and  to  function  on  the 
UNIX  operating  system.  /Assembly  language  particle  movers  have  been 
written  by  L.  E.  Rudslnski  to  maximize  code  performance  on  the  CRAY-1 
computer.  The  code  has  been  run  under  7  different  operating  systems  in 


the  past  decade. 

A  preliminary  input  parameter  list  has  been  constructed  and  run-time 
requirements  for  potentially  useful  problems  are  being  considered. 
Because  it  is  possible  to  overcommit  our  hours  or  computer  time  with  such 
a  code,  we  give  more  immediate  priority  to  improved  modeling  with  the 
pre-existing  1-D  electrostatic  code;  under  our  first  year  effort  we  have 
devoted  most  of  our  small  computational  effort  to  preparing  1-D  code 
modeling. 

But  the  radial  boundedness  of  the  electron  beam  is  expected  to  lead 
to  jetting  of  the  plasma  ions,  and  the  stability  or  overstability  of  the 
radial  ion  motion  cannot  yet  be  predicted.  Radial  oscillations  could 
lead  to  the  emission  of  clearly  defined  bunches  just  as  easily  as  axial 
space-charge  oscillations,  and  the  size  (and  energy  content)  of  the 
bunches  would  likely  be  controlled  by  the  frequency  (etc.)  of  the 
oscillations.  Hence  we  see  2-D  modeling  as  potentially  quite  important 
to  understanding  bunching. 

The  WAVE  code  is  a  relatively  small  code  which  spends  most  of  its 
time  in  just  a  few  small  subroutines.  This  makes  it  relatively  easy  to 
optimize  for  particular  machines  and  to  understand  how  it  works.  The 
code  uses  finite  sized  particles  to  represent  the  electrons  and  ions  and 
solves  Maxwell's  equations  on  a  finite  grid  using  the  sources  of  charge 
and  current  from  the  particles.  At  each  time  step  these  fields  are  used 
to  advance  the  particles  and  the  process  is  repeated. 

The  subroutine  TRANS  obtains  the  particles  from  disk  through  a 
multichannel  triple  buffered  disk  input-output  module.  The  batches  of 
particles  are  brought  in  from  disk  and  sent  to  one  of  three  particle 
moving  routines  depending  on  the  selected  physics  and  to  the  diagnostic 


routine  PDIAG  for  either  movie  or  still  frame  diagnostics  at  the 
appropriate  time.  The  batch  is  then  returned  to  disk  through  the  I/O 
buffer.  If  the  maximum  number  of  particles  will  fit  into  the  I/O  buffer, 
the  particles  are  never  transferred  to  disk..  After  processing  all  the 
particles,  the  electromagnetic  fields  are  advanced  to  the  next  time  level 
using  the  updated  sources  of  charge  density  and  current  density. 
Finally,  at  selected  time  steps  the  DIAGNOS  routine  is  called  to  provide 
graphical  output  of  the  grid  quantities.  A  variety  of  input  controls  of 
these  diagnostics  are  available. 

The  code  requires  a  spatial  mesh  size  small  compared  with  c/w  and 

P 

timesteps  a  small  fraction  of  the  local  plasma  period.  As  a  result,  its 
use  for  macroscopic  or  long-timescale  dynamics  is  likely  to  be 
prohibitively  expensive.  We  expect  it  to  be  useful  as  a  tool  to  simulate 
2-0  behavior  of  space-charge  oscillations  in  the  source  plasma  sheath 
region,  and  perhaps  the  dynamics  of  low-density  background  plasma 
surrounding  a  denser  plasma  blob. 


4.  SOURCE  PLASMA 


In  order  Co  coordinate  closely  with  the  concurrent  experimental 
efforts,  we  analyze  first  the  ionization  processes  in  the  ~  2  cm  gas  puff 
(P  <  .1  torr)  through  which  the  electron  beam  (described  above) 
propagates  and  from  which  it  extracts  ions  to  form  a  partially 
neutralized  propagating  plasmoid  downstream.  As  a  first  approximation  in 
describing  the  electron  and  ion  density  buildup  by  ionization,  we  ignore 
the  effects  of  electric  fields  inside  the  ionizing  gas  puff.  The  gas, 
puffed  transversely  across  the  electron  beam  path  adjacent  to  the  Luce 
diode,  is  exposed  to  beam  electrons  and  cathode  radiation  through  the 
anode  aperture,  approximately  1  cm  in  diameter. 

We  calculate  and  compare  (a)  impact  ionization  by  the  1  MeV 
reflexing  electron  beam,  and  (b)  photoionization  by  cathode  UV 
radiation.  The  beam  is  modeled  as  a  density  n^(t)  given,  for  the  first 
10  ns,  by 

1 0  -A 

n,(t)  =  10  cm  •  (t/10  ns)  . 
b 

Assuming  the  beam  has  spread  radially  to  1  cm  diameter,  this  corresponds 
to  a  peak  one-way  current  of  33  kA  for  a  radially  flat  beam  profile.  For 
H2,  we  expect  this  primary  ionization  to  build  plasma  density  according 
to 
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To  estimate  the  photoionization  rate,  we  suppose  the  cathode  plasma 
to  be  a  blackbody  radiating  as  a  7  eV  plasma.  The  radiation  flux  over  a 
2ir  hemisphere  solid  angle,  from  the  15  mm  radius  cathode  (.0707  cm^), 
with  unit  emissivity,  would  be 


7  T  4 

P(w)  =  1.7  X  lO'Cy-^)  , 


with  a  more  or  less  uniform  radiation  intensity  a  distance  r  (cm)  from 
the  cathode  center, 


T 

yw/cm2)  -  2.7  x  lO6^^-)4  ^ 


The  H2  photoionization  cross-section  rises  quickly  from  the  13.6  eV 

-17  2 

photon-energy  threshold  to  peak  at  cy  =>  10  Mbarns  =  10  cm  (Fig. 

4).  The  integrated  radiation  intensity  above  14  eV  from  a  7  eV  blackbody 

is  about  80%  of  the  total  flux,  so  we  multiply  I,  by  0.8  and  use  the  peak 

<P 

cross-section.  This  gives 
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for  constant  cathode  conditions  during  the  first  10  ns.  (This 

corresponds  to  about  6.5%  ionization  of  the  gas,  increasing  to  ~20%  by 
30  ns  when  the  beam  stops.) 

Thus  we  conclude  that,  even  if  the  emissivity  of  the  cathode  is  as 
low  as  0.1  and  its  temperature  as  low  as  3.5  eV  Instead  of  7  eV, 
photoionization  probably  dominates  over  impact  Ionization  by  the  1  MeV 
beam  electrons.  Tills  is  an  Important  departure  in  the  group's  modeling 


r.r 


of  the  ion-source  plasma  formation. 


The  modeling  of  avalanche  ionization  in  the  electric  fields  of  the 

rising  beam  has  not  yet  been  included.  The  radial  space-charge  electric 

field,  Er,  of  a  .5  cm-radius  beam  with  nb  =  10^cm-^(t/10  ns)  is  of  order 

4  MV/cm  •  (t/10  ns)  at  the  beam  edge,  if  there  is  no  space-charge 

neutralization.  But  in  about  1  ps,  photoionization  can  produce  n^  ~  n^, 

so  that  if  the  Initial  plasma  electrons  are  expelled  by  the  fields, 

charge  neutralization  is  established  very  early  and  no  such  radial  field 

magnitudes  build  up.  Likewise,  inductive  fields  due  to  the  rise  of  net 

current  are  much  reduced  by  the  return  current  carried  by  photoionized 

electrons.  Since  our  model  predicts  n  /n,  >  10^  at  all  times,  these 

P  b  ~ 

return  currents  probably  do  not  greatly  modify  the  plasma  otherwise. 

With  the  preceding  model,  the  ratio  of  beam  density  to  plasma 
density  is  nearly  constant  during  the  beam  rise  phase 
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If  plasma  density  gradients  do  not  interfere,  this  means  a  growth  rate  of 
the  hydrodynamic  electron-electron  two-stream  instability  on  the  order  of 
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10  ns,  the  beam  moves  2.6  cm  In  one  growth  time,  and  so  traverses 


the  plasma  in  about  one  e-folding  time  of  the  unstable  space-charge 
oscillations.)  The  frequency  of  the  unstable  waves  is 
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which  should  be  detectible  by  its  conversion  to  EM  waves  at  this 
frequency  of  order  90  GHz.  The  10  ns  beam-rise  time  represents  3650  or 
so  plasma  oscillation  periods  (their  variation  included),  or  some  67 
growth  times,  and  one  may  thus  expect  a  fully-developed  two-stream 
instability  in  the  source  plasma  unless  gradients  significantly  damp  the 
instability.  (The  plasma  is  only  about  20  wavelengths  thick,  at  the  most 
unstable  wavelength.) 

Such  a  microinstabll ity ,  if  it  exists,  probably  has  two  important 
effects  besides  generation  of  microwave  diagnostic  radiation:  it  can 

make  the  plasma  electrons  effectively  more  collisional  and  less  tied  to 
magnetic  field  lines,  and  it  can  spread  the  energy  distribution  of  the 
beam  electrons.  The  latter  may  already  be  observed,  in  the  large  beam 
temperatures  required  for  consistency  with  the  observed  downstream  beam 
equilibrium,  and  the  former  may  also  be  observable. 

The  photoionization  of  the  10~3  torr  "vacuum"  region  of  the  Maryland 
experiment  may  also  provide  a  space-like  low-density  ambient  plasma 
through  which  the  fast  electrons  and  accompanying  ions  propagate 
downstream.  Here,  according  to  the  propagation  stability  analysis  of 
Sec.  5,  the  growth  rate  in  the  beam  frame  should  be  of  order 
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or  less  if  the  beam  is  now  "hot".  At  3  to  10  cm  from  the  cathode  with 
full  eraissivity,  this  growth  rate  would  vary  from  10^  to  10®s-^  at  t  = 

10  ns,  and  would  have  a  signature  of  microwave  emission  at  and  below 
1.3  GHz  if  present.  This  growth  rate  corresponds  to  .25  to  2.5  m  of  beam 
travel,  and,  since  the  beam  only  traverses  about  1  m  in  the  experiment, 
significant  exponentiat ion  probably  cannot  be  seen. 


5.  STABILITY  OF  PROPAGATING  PLASMOID  -  LINEAR  THEORY 


The  first  step  in  the  theoretical  analysis  of  plasmoid  propagation 
is  to  evaluate  the  stability,  (or  growth  rates,  if  unstable),  of  a  quasi¬ 
neutral  beam  passing  through  a  tenuous  quasi-neutral  plasma.  This  is 
done  by  treating  the  background  plasma  as  a  perturbation  on  the  beam,  in 
the  beam  frame;  i.e.,  reversing  the  usual  roles  of  beam  and  plasma. 

Even  before  a  mostly  force-neutralized  plasmoid  is  'disconnected' 
from  its  source,  streaming  instability  may  affect  its  momentum  scatter, 
if  the  instability  grows  fast  enough  (i.e.,  in  a  time  much  less  than  the 
plasmoid  length  divided  by  the  front  velocity).  But,  for  sufficiently 
tenuous  background,  the  most  important  instability  effects  may  be  seen 
only  on  a  longer  timescale,  once  the  plasmoid  is  propagating  freely. 
This  instability  amplitude  is  enhanced  by  any  modulations  in  density  that 
may  have  come  about  during  plasmoid  formation  and  early  history.  Seen  in 
the  plasmoid  ion  (and  electron)  frame,  the  tenuous  background  streams  by 
axially  at  relativistic  speed  and  causes  the  unstable  growth  of  density 
modulations  (in  both  plasmoid  and  background).  Depending  on  the  absolute 
or  convective  nature  of  the  instability  in  the  plasmoid  frame  (Bers  and 
Ram,  1982),  and  on  the  group  velocity  of  the  perturbations  in  that  frame, 
these  density  modulations  may  pile  up  or  may  propagate  to  the  rear  of  the 
plasmoid  (Kuzelev  and  Rukhadze,  1979).  In  the  latter  case,  if  the  growth 
time  is  not  much  shorter  than  the  convection  time,  the  unstable  wave 
energy  is  shed  into  the  background  behind  the  plasmoid,  and  the  plasmoid 
propagation  may  be  only  weakly  affected  by  the  fields  of  the  instability. 

On  the  other  hand,  if  the  instability  grows  faster  than  it  convects 


to  the  rear  of  the  plasmoid,  nonlinear  processes  probably  set  in  which 


lead  to  quasi-stationarv  density  modulations  (cavitons)  in  the  plasmoid 
frame  (Zakharov,  1972).  It  is  likely  that  these  axial  dependences  of 
density  do  not  break  up  the  plasmoid,  but  they  affect  the  electron 

scattering.  The  issues  (or  two  important  issues)  for  long-range 
propagation  of  such  neutralized  beams  are  (a)  what  is  the  time  required 

for  transverse  momentum  spread,  and  (b)  what  is  the  timescale  on  which 

electrons  (or  ions)  scattered  out  of  the  beam  are  replaced  with 
accelerated  background  particles  (because  this  determines  not  only  an 
energy  loss  rate,  but  also  a  residual  set  of  electric  and  magnetic 

fields).  The  presence  of  instability  microfields  impacts  both  of  these 
issues. 

In  addition  to  the  conventional  streaming  instability  of  space 
charge  bunching,  there  is  a  purely  electromagnetic  f ilamentation 
Instability  (Benford,  1973)  which  can  strongly  enhance  the  rate  of 
transfer  from  parallel  to  transverse  momentum.  For  very  energetic  beams 
in  axial  magnetic  guide  fields  this  instability  can  grow  even  faster  than 
the  electron-electron  two-stream  (Ottinger  and  Guillory,  1979),  although 
for  moderate-gamma  beams  It  is  more  slowly  growing.  Even  if  this  mode  is 
considerably  slower  than  two-stream  modulation  growth,  it  may  affect  the 
transverse  beam  spread  more,  because  of  the  transverse  nature  of  the 
unstable  fields  and  the  different  wave  convection  rate  along  the  beam. 
We  have  not  treated  this  f ilamentat ion  mode  in  our  first-year  effort,  but 
have  turned  our  attention  to  the  electron-electron  two-stream  interaction 
discussed  above.  The  ID  ES  model  is  shown  in  Fig.  5  and  the  main  results 
from  this  Initial  analytical  investigation  are  tabulated. 


6.  STABILITY  OF  PROPAGATING  PLASMOID  -  NONLINEAR  THEORY 


Computer  simulation  of  the  nonlinear  stage  of  electron-electron  two- 
stream  instability  in  a  bounded  system  with  the  beam  exit  boundary  "free" 
have  been  ongoing  in  the  University's  Space  Sciences  Laboratory  under  NSF 
funding  and  were  made  available  without  charge  to  this  contract.  We 
report  here  some  of  the  results  of  those  simulations  by  H.  Rowland. 

A  1-D  electrostatic  particle-in-cell  code  was  run  with  a  warm 
tenuous  'beam'  (b)  entering  from  the  left  of  a  region  of  initially- 
uniform,  equally  warm  plasma  (p).  In  the  run  shown  here,  loss  of  plasma 
particles  to  the  boundaries  was  compensated  by  injection  of  Maxwellian 

plasma  from  the  right-hand  boundary,  and  injection  of  'beam'  plus 

Maxwellian  distribution  from  the  left.  The  right-hand  boundary  potential 
was  allowed  to  float.  The  ratio  of  densities  n^/np  for  the  run  shown  was 
0.02. 

The  code  was  run  both  with  fixed  (infinitely  massive)  ions  and  with 

mi/me  =  100  and  T^  =  Tg.  By  about  600  a  relatively  stationary 

soliton  of  width  30-40  Ap  was  observed  at  x  =  180  from  the  injection 

side.  Its  amplitude  was  about  3  1/2  times  larger  with  mobile  ions  than 

with  fixed  ions.  Both  cases  are  shown  in  Figs.  6a  and  6b  at 

t  =  650  a)  ,  Vertical  units  are  arbitrary  (i.e.,  not  normalized)  and  the 
P 

2 

vertical  scale  is  proportional  to  E  /8w . 

With  m^/me  =  100,  the  field  tended  to  accelerate  some  plasma 
electrons  antiparallel  to  the  beam  direction,  as  shown  by  the  logarithmic 
plot  of  the  parallel  velocity  distribution  (b+p)  at  650  w  ^  in  Fig.  6b. 
The  flux  of  suprathe  rma  1  electrons  at  the  boundary  tended  to  have  an 
erratic  spiky  time  hlstorv.  These  results  were  presented  by  Rowland  at 


the  AGU  meeting  in  Baltimore  in  May  1986  and  further  studies  were 
described  at  the  November  1986  APS-DPP  meeting. 

We  have  suggested  further  runs  of  this  code  with  a  cold  beam  and 
with  modified  plasma  particle  boundary  conditions. 

2 

Calculations  of  the  peak  electrostatic  oscillation  energy  E  /8u  . 

max  ’ 

normalized  to  the  background  streaming  energy  in  the  plasmoid  frame 
2 

(Y-l)n^mc  ,  for  Y  <  10  (e.g.,  Matsiborko  et  al. ,  1971)  have  shown  that 


W  =  E^ax/[8TTnbmc2(Y-l)  ]  ~  (Y  + 1 )  (nfe/n  )  1  ^  3/ Ax 


when  Y(n, /n  )  <<  1  and  when  saturation  can  occur  inside  the  plasmoid. 

b  p 

For  a  background  density  n^  =  lO^cm-^,  Y  =  2,  and  n0  =  101(^cm-^  as 

-92  -23 

an  example,  one  has  W  =  10  ,  E  / 8mr  =  0.9  *  10  erg/cm  and  E  ~  190  V/cm 

at  saturation.  The  ponderomotive  pressure  in  the  field  peak  corresponds 
to  only  7  *  10  ^torr,  compared  with  a  thermal  pressure  in  the  plasmoid 
frame  of 


P^Ctorr)  ~  1.2  x  10  (  1QP  -3)  ( io'o'PeV^  * 

10  cm 


Thus  we  do  not  expect  hydrodynamic  disruption  of  the  plasmoid  by 
'solitons'  in  such  a  low-density-background  case. 

But  for  n^  =  lO^cm-^  (as  could  be  produced  by  photoionization  of 
10--*  torr  laboratory  vacuum  in  a  drift  tube),  the  ponderomotive  pressure 
at  saturation  is  about  500  times  larger  and  is  comparable  to  the  thermal 
pressure  in  the  plasmoid,  leading  on  long-enough  timescales  to  local 
expulsion,  i.e.,  cavitation,  of  the  plasma.  The  required  timescales, 
however,  would  not  be  seen  in  tvpical  drift  tubes;  for  y  =  2,  saturation 


(~  10  e-folds)  f  the  instability  requires  about  6.5  meters  at 
10®cm~^,  and  hydrodynamic  cavitation  requires  about  20  meters. 
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INTENSE  BEAN  PROPAGATION  PROPERTIES  IN  MACNETI2ED  AND  LOCALIZED  ION  SOURCE  CONFIGURATIONS 


W.  W.  Descler,  R.  J.  Faehl,  P.  G.  O'Shea,  M.  Reiser,  Z.  Segalov,  C.  D.  Strlffler,  and  X.  Zhang 
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Abstract 

During  the  past  few  years  we  have  performed 
both  experimental  and  theoretical  studies  on  the 
propagation  of  intense  electron  beams  Into  an 
evacuated  drift  tube.  In  one  set  of  studies  the 
radial  force  balance  is  achieved  by  an  applied 
magnetic  field  and  in  a  second  set  of  studies  the 
radial  force  balance  is  achieved  by  ions  that  have 
been  created  in  a  localized  region  of  the  tube  and 
collectively  accelerated. 

Introduction 

In  recent  years,  the  propagation  of  intense 
relativistic  electron  beams  (IREB's)  has  been  the 
subject  of  many  studies,  in  part  due  to  their 
various  applications  in  such  diverse  fields  as 
high  power  coherent  radiation  source  development, 
collective  ion  acceleration,  and  plasma  heating 
and  confinement.  A  summary  of  the  field  prior  to 
1982  is  presented  in  the  book  by  R.  B.  Miller.* 

A  schematic  of  the  experimental  system  in  our 
studies  Is  shown  in  Fig.  1.  The  IREB  is  field 
emitted  from  a  cathode  1.2  cm  upstream  of  a 
stainless  steel  anode.  The  nominal  properties  of 
the  electron  beam  Injected  into  the  vacuum  chamber 
are  1  MeV,  25  kA,  30  ns  FWHM.  Various  diagnostic 
equipment  is  employed  in  the  experiments,  one  of 
which  is  the  current  collector/caloriraeter  shown 
in  the  figure.  In  the  first  set  of  studies 
discussed,  a  magnetic  field  immerses  the  entire 
system  and  in  the  second  set  of  studies,  a  puff 
valve  supplying  a  localized  gas  cloud  is  placed 
immediately  after  the  anode.  No  magnetic  field  is 
applied  in  the  latter  studies. 


CURRENT  COLLECTOR  /  CALORIMETER 


FIG.  1.  Experimental  configuration. 


Magnetically  Confined  Beam  Propagation 


For  the  magnetically  confined  studies,  the 
propagated  current  is  determined  as  a  function  of 
the  uniform  applied  magnetic  field  (0-2  T)  for 
various  drift  tube  diameters,  cathode  geometries, 
and  anode  aperture  sizes. 

Experiments 

The  experimental  configuration  is  shown  in 
Fig.  I  with  the  absence  of  the  puff  valve.  A 
solid  cathode  and  anode  apertures  of  1.2  and 
2.6  cm  are  used  In  combination  with  stainless 
steel  drift  tubes  of  3.8,  9.8,  and  14.8  cm  Inner 
diameter.  The  beam  current  propagating  to  the  end 
of  the  drift  tube  Is  measured  using  a  low 
Impedance  current  collector.  The  injected  beam 

current  la  measured  by  placing  the  collector  2  cm 
downstream  of  the  anode  plane  and  is  found  to  be 
27  *  1  kA  Independent  of  applied  magnetic 


field. 

10  -10 


^Tie  vacuum 
Torr. 


maintained 


la 


in  the  range 


The  peak  electron  beam  current  measured  at 
the  downstream  end  of  the  70  cm  long  drift  tube  as 
a  function  of  applied  magnetic  field  B  is  shown  in 
Fig.  2  for  the  various  geometries  considered.  The 
following  Important  features  are  apparent  from  the 
measurements: 


1.  In  all  tubes  there  appear  to  be  three 
distinct  regimes  of  beam  propagation.  These 
distinctions  are  most  obvious  In  the  14,8  and 
9.8  cm  diameter  tube  results. 

(a)  Low  Magnetic  Field.  Propagating  current 
Initially  rises  a6  B  is  Increased  from 
zero.  Here  we  have  a  poorly  confined  beam,  a 
significant  fraction  of  which  may  hit  the 
tube  wall.  As  B  Increases,  the  beam 
confinement  Improves  and  more  current 
propagates . 

(b)  Intermediate  Magnetic  Field.  Beam  current 
decreases  rapidly  with  increasing  B  with 
space-charge  effects  playing  an  increasingly 
significant  role  In  this  region. 

(c)  High  Magnetic  Field.  A  plateau  region  where 
the  change  in  propagating  current  Is  small  as 
B  increases. 

2.  As  tube  diameter  decreases,  the  region  of 
maximum  current  propagation  shifts  to  higher 
magnetic  fields,  and  the  propagating  current 
in  the  high  field  plateau  region  increases. 

3.  For  the  14.8  and  9.8  cm  diameter  tubes, 

changing  the  anode  aperture  size  from  1.2  to 
2.6  cm  diameter  did  not  greatly  alter  the 
propagating  current.  However,  we  note 
significant  difference  In  the  current 
characteristics  for  the  high  field  region  for 
the  3.8  cm  tube.  This  is  due  to  the  fact 

that  the  anode  diameter  is  a  significant 
fraction  of  the  tube  diameter.  The  smaller 
anode  hole  case  appears  to  exhibit  similar 
features  in  the  region  0.8  <  B  <  1.1  T  as  the 
large  tube3  do  in  the  range~0.3~<  B  <  0.5  T.  * 
We  have  evidence  that  the  current  propagating 
In  this  tube  for  the  2.6  anode  begins  to 
decrease  at  B  ■  1.2  T  and  above. 

The  fact  that  the  smallest  tube  can  propagate 
as  much  current  at  high  B  as  the  large  onea  can  at 
considerably  lower  B  removes  from  consideration 
the  possibility  that  some  form  of  magnetic 
mirroring,  resulting  from  nonuniformity  of  the 
applied  field  causes  the  observed  reduction  of 
current  in  the  large  tubes  at  high  magnetic  fl^ld. 

Model  of  Solid  Beam  Equilibrium 

A  steady-state  model  of  an  IREB  in  an  axla- 
aymmetric  system  can  be  constructed  by  conserving 
single  particle  energy  and  canonical  angular 
momentum,  applying  continuity  of  current,  and 
maintaining  force  balance.  The  beam  model  and  Its 
various  parameters  are  shown  In  Fig.  3.  We  assume 
that  a  solid,  uniform  density  electron  beam  of 
radius  R  la  Injected  into  a  long  grounded 
cylindrical  drift  tube  of  radius  R  .  The  injected 
beam  is  monoenergetlc  with  energy 

iic2(r  -  1)  -  «v  . 

o  o 


After  the  bean  pt.se.  through  the  anode.  It 
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expanda  ad labac Ically  to  a  radius  where  a 

laminar  flow  equilibrium  is  set  up  far  from  the 
end  walls.  The  entire  system  is  Immersed  in  a 
uniform  magnetic  field  B^ .  The  downstream  beam 
properties  are  the  charge  density  p(r)  and  the 
azimuthal  and  axial  velocities  V  (r)  and  V  (r). 
The  fields  generated  by  the  beam  are  confined 
inside  the  drift  tube,  and  include  the  radial 
electric  field  E  (r)  and  the  azimuthal  and  axial 

magnetic  fields  I?  ,  and  B  . 
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FIG.  2.  Peak  propagating  current  measured  at  the 
current  collector  for  various  drift  tubes  and 
anode  apertures  using  a  solid  cathode. 


To  further  simplify  the  analysis,  the 
equilibrium  beam  density  Is  assumed  constant 
across  the  beam  p(r)  -  p^,  allowing  the  electric 
potential  +(r)  and  the  electric  field  to  be  easily 
calculated.  Conservation  of  single  particle 
energy  gives 


l(r)  -  y0  +  -Sj  *(r)  -  [1  -  B?(r)  -  g]<r)j",/2  , 

me 

Cl) 

where 

pbRh  r2  R, 

♦(<•)  -  [i  -*1  *  i  im~)]  .  <2) 

o  Rb  % 


In  a  laminar  flow  beam  of  uniform  density, 
canonical  angular  momentum  conservation  can  be 
written  as 


eB  ?t  R  2 

• - 2 —  (a-)  "  rl<aT<f>VACr)  -  eAA(r)]  , 


V 


(3) 


where 


A  ( r )  -  rB 

i  ♦ 


./ 2 


A  (r) 


, . ,  is  the  total 

azimuthal*  compon«A?t  of  8*Che  magnetic  vector 
potential.  Solving  Eq  .  (3)  for  V  (r)  and 

substituting  It  into  J  (r)  -  p  V  (r),  le  obtain 
from  Maxwell's  Equations  an*  inhomogeneous 
differential  equation  for  A  ,  This  equation  can 
be  aolvcd,  from  which  R  is'^computed,  then  V  (r) 
Is  calculated  by  Eq .  (Ujt  V  (r)  via  Eq.  (l),*and 
B  from  Ampere's  Law,  2 
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FIG.  3.  Schematic  of  the  beam  model. 


At  this  point,  the  functional  forms  of  all 
the  equilibrium  beam  properties  have  been  found  in 
terms  of  R  ,  R  ,  R  ,  y  ,  B  ,  and  the  beam  density 
0  •  In  addition, wwe  ^equiTe  each  volume  element 
of°the  beam  to  be  in  radial  force  balance,  which 
can  be  written  as 


PY(r)V^(r) 


-  e(Esr<r>  +  Vr>[BAZ  +  Bsz(r)] 


-  Vr)Vr)1  * 


(4) 


To  solve  this  equation,  R  ,  R  ,  R  ,  y  ,  and  B 
are  held  fixed  while  p^  lsa  varied  Ji\tll0the  "bel? 
possible  equilibrium0  is  found.  After  calculating 
the  beam  density  for  an  equilibrium,  the  axial 
beam  current  19  calculated  from 


I  -  2„p  /  V  (r)  rdr  ,  (5) 

6  0  * 

which  equals  the  Injected  current. 

From  the  above  model,  we  can  generate  the 
maximum  propagated  current  as  a  function  of 
applied  magnetic  field  B^.  A  specific  example  is 
shown  In  Fig.  A.  In  the  limit  of  large  magnetic 
field,  R.  ♦  R  ,  the  peak  beam  current  agrees  well 
with  the  fiogdankevlch-Rukhadze  limiting  current2 


BR 


17(y2/J  -  I)3/2 

l  +  2  in( R  /R  ) 

w  b 


IkA) 


(6) 


which  is  plotted  for  reference.  As  B  decreases, 
the  beam  expands  (1L>  R  )  and  the  pe5?  current  Is 
limited  by  space  charge  Effects,  A  rapid  Increase 
occurs  In  the  region,  0.4-0. 6  T,  where  the  current 
i9  limited  by  total  field  energy.  The  beam  fills 
the  tube  below  the  peak,  -  0.3  T,  where  the 
current  is  limited  by  the  size  of  the  tube. 

In  comparison,  the  model  predicts  all  the 
experimental  trends.  The  main  disagreement  is  in 
the  magnitude  of  the  applied  magnetic  field.  Here 
the  theory  is  off  by  approximately  a  factor  of 
two.  It  is  expected  that  cyclotron  motions,  other 
nonlamlnar  effects,  and  nonadlabatlc  expansions 
will  account  for  some  of  this  discrepancy.  A 
detailed  description  of  this  system  is  presented 
in  Ref.  3. 

Localized  Ton  Source  Ream  Propagation 

For  the  localized  ion  source  studies,  wc 
measure  the  propagated  current  as  a  function  of 
the  Ion  source  properties.  The  ion  source  Is  from 
a  localized  gas  cloud  In  the  immediate  vicinity  of 
the  injected  electron  beam  (puff  valve).  A 
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virtual  cathode  forms  In  the  Ion  source  region, 
where  the  strong  electric  field  (100-500  MV/ra)  la 
able  to  accelerate  the  ions  downstream  and  allow 
the  formation  of  a  propagating  beam.'1* 
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baz.  axial  applied  MAGNETIC  FIELD  (T) 

FIG.  4.  Maximum  beam  current  versus  applied 
magnetic  field  for  two  different  anode  sizes. 

R  “  1.9  cm,  y  “  3.5. 
w  o 

Experiments 

The  experimental  configuration  is  shown  in 
Fig.  1.  An  IREB  is  emitted  from  a  6  nan  diameter 
tungsten  cathode.  A  25  mm  aperture  In  the  anode 
allows  almost  all  of  the  beam  electrons  to  pass 
through  into  the  downstream  drift  region.  The 
drift  tube  diameter  Is  15  cm.  Ions  to  be 
accelerated  are  produced  by  electron  impact  and 
lon-lon  avalanche  ionization^*  of  the  localized  gas 
cloud.  The  firing  of  the  fast  gas  puff  valve  is 
adjusted  to  ensure  that  the  gas  cloud  Is  localized 
within  2-3  cm  of  the  anode  plane  at  the  time  of 
electron  beam  injection. 

The  diagnostics  used  are  a  number  of  current 
probes  to  measure  both  the  axially  propagating 
current  and  the  radial  current  at  the  drift  tube 
walls.  In  addition,  a  7.5  cm  diameter  combination 
Faraday  cup/calorlmeter  measures  the  Integrated 
net  axial  beam  current  and  energy. 

Measurements  of  the  net  current  propagated  to 
a  given  axial  position  using  the  Faraday  cup  are 
displayed  In  Fig.  5(a).  Shown  In  Fig.  5(b)  Is  the 
energy  deposited  In  the  calorimeter  and  in 
Fig.  5(c)  the  counts  from  a  silver  activation 
neutron  detector  as  a  function  of  gas  cloud 
pressure.  When  the  peak  pressure  of  the  localized 
gas  cloud  Is  adjusted  to  an  optimum  value  (about 
35  mTorr),  beam  current  many  times  the  B-R  space 
charge  limiting  value,  Eq.  (o),  can  be  propagated 
down  the  drift  tube.  In  the  absence  of  the  Ion 
source,  currents  well  below  the  space  charge  limit 
are  observed.  At  the  optimum  value,  a  large 
fraction  of  the  Injected  current  propagates  to  the 
axial  location  of  the  Faraday  cup.  The  range  In 
pressure  about  the  optlmlum  narrows  as  the  cup 
location  Is  moved  farther  downstream.  This 
appears  to  be  limited  only  by  the  pulse  duration 
of  the  Injected  electron  beam  or  by  the  depletion 
of  the  source  of  Ions  at  the  Injection  point. 

Using  a  radial  array  of  axial  current  probes , 
measurements  have  been  made  of  the  propagating 
current  profile.  Several  features  are  worthy  of 
mention:  (l)  the  current  profile  Is  virtually 
flat  and  the  measured  current  levels  are  very  low 
when  no  localized  ion  source  is  present,  (2)  the 
best  radial  confinement  of  the  beam  electrons  is 
observed  at  a  peak  gas  cloud  pressure  of  about 
35  mTorr  (corresponding  to  the  optimum  pressure 
for  peak  net  current  propagation  discussed  above). 


and  (3)  a  hollowing  out  of  the  beam  la  observed 
when  the  gas  cloud  pressure  exceeds  this  optlmlum 
value. 


AP  (RELATIVE  LTaTS) 

FIG.  5.  Results  of  (a)  current  collector, 
(b)  calorimeter,  and  (c)  neutron  detector 
measurements  at  two  axial  locations  as  a  function 

of  gas-cloud  peak  pressure  (p  ). 

max 

An  axial  array  of  wall  probes  Installed  to 
measure  the  radial  current  flowing  to  the  walls 
has  resulted  In  considerable  Insight.  Several 
features  are  again  worthy  of  note:  (l)  with  no 

localized  ion  source  at  the  injection  point,  the 
wall  currents  are  quite  high  at  the  Injection  end 
and  fall  off  to  very  low  values  as  the  probe  axial 
position  Is  Increased;  (2)  for  the  probe  closest 
to  the  Injection  point  the  observed  current  Is 
lowest  at  a  peak  pressure  of  about  30  mTorr,  the 
same  pressure  at  which  the  high  downstream  vail 
currents  are  observed;  and  (3)  gas  cloud  pressures 
above  the  optimum  result  In  higher  wall  currents 
at  the  Injection  end.  These  as  well  as  a  more 
extensive  presentation  of  the  experimental  studies 
appear  In  Refs.  5  and  6. 

Model  of  Electron-Ion  Beam  Equilibrium 

A  steady-state  electron-ion  beam  model  la 
formulated  that  relates  the  downstream  equilibrium 
parameters  to  the  Initial  Injection  properties  of 
the  system.  Because  there  Is  no  applied  magnetic 
field,  the  propagation  depends  on  the  presence  of 
Ions  which  come  from  a  localized  source  near  the 
entrance  end  of  the  Injected  beam.  The  electron 
beam  Is  assumed  to  have  a  transverse  temperature 
T  In  addition  to  Its  axial  velocity  and  the  Iona 
ha$e  only  an  axial  velocity.  A  schematic  of  the 
model  Is  shown  Lp  Fig.  6.  An  electron  beam  of 
kinetic  energy  «c  (y  -  1)  -  eV  and  current  I  la 
injected  into  a  long  cylindrical  drift  tube0  of 
radius  R  .  The  ion  source  parameters  ire  a 
format ionwrate  S  (C/sec)  in  a  potential  (V). 
As  shown,  the  iyntea  evolves  into  a  downstream 
equilibrium  that  contains  both  ions  and  electrons 
with  axial  velocity  of  v  -  8  c  a'd  *  “Be, 
respectively,  in  a  poteiuial  v.  The  Vy»  teatein 
general  possesses  a  net  charge  and  net  current. 

The  downstream  equilibrium  of  the  Ions  and 
•lectroni  requires  radial  force  balance  which 


results  la  a  radial  density  profile  for  the 
electrons  and  Ions  of 


"e(r>  -  - - 77  -  2.2  -  T~&  V  ■  nl(r)  (?) 

(1  +  (r/a)  ]  zl  ze 


8c  kT  1  -  B  . 

n  -  -2-7 - ~ ■  .  (8) 

e  a  (Bze  -  Bzl)2 

This  profile  is  characterized  by  a  peak  density  a 
and  an  effective  radius  "a."  Continuity  o? 
current  for  each  species  gives 
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Conservation  of  energy  gives 
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where  fl2  «  kTe/">ec2.  The  on-axla  downstream 
potential  Is  ■ 

2kT  B  R  , 
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ze  zi 

and  the  net  current  can  be  calculated  from 

1  *  s<  -  I  . 

net  1  o 
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Electron  Ion  Equilibrium 

Beam  Source 


V  I  ■ 


FIC.  6,  Schematic  of  the  model* 


The  above  model  is  analyzed  for  the  case 

I  -  25  kA,  eV  -  l  MV,  and  R  /a  -  10.  An 
electron  beam  temperature  is  chosen  and  then  the 

above  equations  are  solved  for  B  ,6  ,  V,  , 

and  5^ .  These  results  are  shownZCin  ffg.  7  whc?e 
we  have  plotted  normalized  net  current  versus 
electron  beam  temperature.  The  regime  of  Interest 
to  ua  Is  where  I  / 1  is  near  -  1.  When  it  Is 

exactly  -  1 ,  T  -nCl**  *°  240  keV  and  B  *  0.  That 

Is,  the  system  ela  charge  neutralized  and 
represents  the  traditional  Bennett  equilibrium. 
The  regime  In  temperatures  Ju9t  below  TJ, 

200  kcV  £  <  TJ  "  240  keV,  gives  rise  to  • 


system  that  is  almost  charge  neutralized 
[f  H  n^/n  >  0.96),  the  ion  source  rate  la  small 
[S./I  <  8.10],  the  ion  potential  at  the  source 

is  smail"*  (Vt  /Vo  <  2.0  going  to  zero  at  T* )  ,  and 
the  ion  axla¥  velocity  Bzi  <  0.07.  These  Results 
are  in  good  agreement  with  the  experimental 
results  and  the  source  rate  is  in  excellent 
agreement  with  ionization  rates  for  the  localized 
neutral  gas  source. 
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FIG.  7.  Normalized  net  current  versus  electron 
beam  temperature.  The  system  parameters  are; 
I  -  25  kA,  eV  -  1  MV,  and  R  /a  -  10. 


Summary 

We  have  performed  both  experimental  and 
theoretical  studies  on  the  beam  propagation 
characteristics  of  two  systems,  one  possessing  a 
magnetic  guide  force  and  the  other  an  ion  guide 
force.  Each  system  contains  a  regime  where  a 
substantial  portion  of  the  injected  electron 
current  is  able  to  propagate  to  distances  far  from 
the  injection  end.  This  current  ie  much  larger 
than  the  B-R  space  charge  limited  current.2  Our 
theoretical  models  agree  well  with  the 
experimental  observations  not  only  in  a 
qualitative  sense  but  for  many  system  variations. 
In  a  quantitative  sense. 

Acknowledgments 

This  work  is  supported  by  the  U.S.  Air  Force 
Office  of  Scientific  Research  and  the  U.S. 
Department  of  Energy. 

References 

1.  R.  B.  Miller,  An  Introduction  to  the  Phvslca 

of  Intense  Charged  Particle  Beans,  (Plenum, 
New  York ,  1982).  “ 

2.  L.  S.  Bogdankevlch  and  A.  A.  Rukhsdzc,  Usp. 

Fir.  Nauk.  101,  609  (1971)  (Sov.  Phya.  Uap. 
_14_,  163  (  19  71TT. 

3.  P.  C.  O'Shea,  D,  Welsh,  W.  W,  De6tler,  and  C. 
D.  Strlffler,  J.  Appl.  Phya.  _55,  3934  (1984). 

4.  H.  Dsntaker,  L.  E.  Floyd,  and  C.  D. 

Strlffler,  Ball.  An.  Phya.  Soc.  29,  1333 

(1984). 

5.  w.  W.  Dcstler,  P.  C.  O'Shea,  and  M.  Reiser, 

Phya.  Rev.  Lett.  32,  1978  (1984);  Phya. 

Flulda  22_,  1897  (  198»T. 

6.  W,  W.  Dustier,  P.  G.  O'Shea,  and  Z.  Scgalov, 
IEEE  Trans.  Nucl.  Scl.  NS- 12 ,  3481  (1985). 


Analytic  Model  of  Intense  Beam  Propagation  through  a 
Localized  Source  of  Ions  into  Vacuum.*  X.  ZHANG  and  C*  D. 
STRIFFLER,  University  of  Maryland. -~A  model  of  the  propagation 

OF  AN  INTENSE  RELATIVISTIC  ELECTRON  BEAM  INTO  VACUUM  AFTER 


PASSAGE  THROUGH  A 
PRESENTED* ^  IT 


LOCALIZED  PLASMA 


BEEN  PREVIOUSLY 


PRESENTED*  J  IT  INVOLVES  A  STEADY“STATE  SYSTEM  COMPOSED  OF 
ELECTRONS  WITH  FINITE  TRANSVERSE  TEMPERATURE  AND  COLD  CO-MOVING 

ions*  Downstream  radial  force  balance,  conservation  of  energy, 

AND  CONTINUITY  OF  CURRENT  LED  TO  RESULTS  FROM  THIS  MODEL  THAT 

O 

AGREE  WELL  WITH  EXPERIMENTAL  OBSERVATIONS.  In  THIS  PAPER,  WE 
PRESENT  RESULTS  OF  THE  MODEL  THAT  INCLUDE  REALISTIC  PULSE 
SHAPES,  FINITE  TRANSIENT  TIMES,  AND  LIMITATIONS  ON  ION 
FORMATION  RATE  DUE  TO  GAS  CLOUD  PARAMETERS*  In  ADDITION,  NEW 
VERSIONS  OF  THE  MODEL  ARE  PRESENTED  AND  ANALYZED* 

*THIS  WORK  IS  SUPPORTED  BY  AFOSR  AND  DOE* 

^Z*  Segalov,  et  al*.  Bull*  Am*  Phys*  Soc*  30 ,  1505  (1985). 

^W*  W.  Destler,  et  al*.  Sixth  Int*  Conf*  on  High-Power 
Particle  Beams,  June  9-12,  1986,  Kobe,  Japan* 


ZY.V.'V 


REFERENCES 


1.  R.  B.  Miller,  "An  Introduction  to  the  Physics  of  Intense 
Charged  Particle  Beam”,  (Plenum,  New  York,  1982) 

2.  R.  J.  Faehl,  X-S,  LAXL,  memorandum,  Nov.  21,  (1983) 

3.  W.  W.  Dcstlcr,  ct  al.  Pliys.  Rev.  Lett.  52.(1984) 

4.  H.  Dantskcr,  ct  al.  Bull.  APS  J29.  1353  (1981) 

5.  Z.  Segalov,  ct  al.,Bull.  /VPS  30,  1505  (1985) 

0.  W.  W.  Dcstlcr,  ct  al..  Sixth  Int.  Conf.  on  High-Power  Partl- 
cal  Beams,  June  9-12,  1980.  Kobe,  Japan. 


I.  Introduction 


II.  Analysis  a)  Radial  Force  Balance 

b)  Continuity  of  Current 

c)  Conseryatlon  of  Energy 

III.  Results  a)  T,  =o 

b)  T,  7^0 


IV.  Application  to  Realistic  Pulse  Shape 


V.  Conclusions 


A  Steady-State  Beam  Propagation  Model 
for  an  Electron  Beam  Injected 


through  a  Localized  Plasma  System 


Injected 

Electron 

Beam 


Localized 
Ion  Source 


Downstream  Equilibrium 
Steady  State  3/3z  =  0 


0 - —  ze 

0  —  zi 


R 

w 


■  iUi  lift 


ANALYSIS 


1.  Downstream  Equilibrium - Radial  Force  Ba 

ance 


Equations  : 

dn; (r ) 
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7’  ,  T - transverse  temperature  of  ion  :ind  electron. 


Solutions  : 
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2.  Continuity  of  Current 


Equations  : 


ions:  St  —  f  ezi  ni  (r  )vzi  ds  =  2ircc  f3:i  J  z-n-(r)r  dr 
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using  the  results  of  force  balance  ,  obtain 
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3.  Conservation  of  Energy 


Equations  : 
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Solution  Procedure 


Given  I0,  VQ, — —  ,Te  and  Tt ,  can  solve  for  (3ze  ,j3:i  ,St  and  Vu 
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Single  Electron  Motion  in  Electric-Magnetic  Field  Generated  by  Fluid 

Motion  equations: 

dp  r  •  • 
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Let  r  =  rQ  f  rr  If  one  linearizes  these  equations  about  rQ 
and  assumes  r1«a,  he  arrives  at  a  solution  with  a  single  fre¬ 
quency  of  oscillation  which  depends  on  r Q.  However,  the  ampli¬ 
tude  of  the  displacement,  r  lmax  Is  >  a .  This  means  the  magnetic 
field  Is  too  weak  to  confine  the  electron  motion  Inside  the  beam 
radius  If  only  taken  to  linear  order.  Therefore,  we  must  keep  high 
order  terms. 
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Fractional  Charge  Neutralization  - Relative  Potential  Depression 


SPECIAL  CASES  Tt  =  c 


1.  CHARGE  and  CURRENT  NEUTRALIZED  LIMIT 
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Realistic  Pulse  Shape  Results 


From  simulation  •  results  of  gas  cloud  ionization  [Dantsker  etal , 
Bull.  APS  29,  1353(1084)1: 


— - -<0.1  and 
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Using  .  these  results  and  (3  ■  >0,  can  determine  the  realistic  pulse  shape 


Conclusions 


1.  A  self-conslstant  downstream  equilibrium  exists  between 
electrons  and  Ions  with  no  axial  magnetic  held.  The  electrons  are 
confined  by  the  self  magnetic  field  and  the  Ions  by  the  self  electric 
field. 

2.  Ti  =0  A  temperature  T t  =  T*  (0zi  =0)  exists  where 
/  =1  and  Inet=I0  (charge  neutral-Bennett  system).  A  parameter 
regime  In  Te  exists  where  the  system  is  nearly  charge  neutralized 
and  the  net  current  Is  nearly  the  Injected  electron  beam  current 
for  a  given  I0  VQ  .  This  regime  Is  Te  <  T*.  By  Increasing  the 
diode  voltage  and  the  Injected  current  T*  will  move  to  higher 
values. 

3.  T,-7^0  The  effect  of  Ion  temperature  T,-  makes 
Te*  (Pzi~ 0)  move  to  lower  values.  Furthermore,  there  Is  no  charge 
neutral  and/or  current  neutral  points  In  regime  Te<Te. 

4.  The  realistic  transmitted  pulse  shapes  change  with  elec¬ 
tron  temperature  and  Ion  temperature.  Pulses  become  narrower 
but  allow  more  net  current  when  temperature  Increases  up  to 
T*(Im  )  —  no  current  can  be  transmitted  ami  Injected  current  Is 


totally  rcllectcd. 


